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Abstract 

We study singularity formation in spherically symmetric solitons 
of the (4+1) dimensional Yang Mills model and the charge two sector 
of the (2+1) dimensional S 2 sigma model, also known as CP 1 wave 
maps, in the adiabatic limit. These two models are very similar. Stud- 
ies are performed numerically on radially symmetric solutions using 
an iterative finite differencing scheme. Predictions for the evolution 
toward a singularity are made from an effective Lagrangian and con- 
firmed numerically. In both models a characterization of the shape of 
a time slice f(r,T) with T fixed is provided, and ultimately yields an 
new approximate solutions to the differential equations that becomes 
exact in the adiabatic limit. 
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1 Introduction 



In this paper we study two hyperbolic partial differential equations that 
develop singularities in finite time. 

The Yang Mills model displays fast blow up, and the S 2 sigma model 
displays both slow blow up and fast blow up. In slow blow up, all relevant 
speeds go to zero as the singularity is approached. In fast blow up the relevant 
speeds do not go to zero as the singularity is approached . The charge 1 sector 
of the S 2 sigma model, investigated in || and |3J, exhibits logarithmic slow 
blow up. The charge 2 sector of the S 2 sigma model and the similar charge 
1 sector of the Yang Mills (4+l)-dimensional model both exhibit fast blow 
up. 

The Yang Mills Lagrangian in 4 dimensons is a generalization of Maxwell's 
equations in a vacuum, and is discussed at length in |IJ. We can regard this 
problem as being that of a motion of a particle, where we wish our particles 
to have certain internal and external symmmetries, which give rise to the 
various gemoetrical objects in the problem. The states of our particles are 
given by gauge potentials or connections, denoted A, on M 4 , and we identify 
M. 4 as H, the quarternions, x = xx+xzi+xzj +x±k. The gauge potentials have 
values in the Lie albegra of SU(2) which can be viewed as purely imaginary 
quarternions, Im(H). The curvature = diAj — djAi + [A^ Aj], where 
[A4, Aj] = AiAj — AjAi is the bracket in the Lie Algebra, gives rise to the 
potential V(A) = (F, F) which is a nonlinear function of A. The static action 
is: 

L = ±j (Fij, Fij)d 4 x. (1) 

The local minima of (0) are the instantons on 4 dimensional space. These 
correspond to solutions of Maxwell's equations in the vacuum. We now 
consider the wave equation generated by this potential with Lagrangian: 

L = \J (d t A, d t A { ) - l -(F l31 F, t] )d A x. (2) 

Via the calculus of variations, the evolution equation for this (4+l)-dimensional 
model is 

d\A i = -VjF ir (3) 

Theoretical work on the validity of the geodesic or adiabatic limit ap- 
proximation for the monopole soltions to the Yang-Mills-Higgs theory on 
Minkowski space is presented in ||. 
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The two-dimensional S 2 sigma model, especially the charge 1 sector, has 
been studied extensively over the past few years in 0, |J, @,|7j], @, JT0[ 



It is a good toy model for studying two-dimensional analogues of elementary 
particles in the framework of classical field theory. Elementary particles 
are thereby described by classical extended solutions of this model, called 
solitons. This model is extended to (2+1) dimensions. The previous solitons 
are static or time-independent solutions, and the dynamics of these solitons 
are studied. The charge 2 sector of this model is held by general wisdom to 
be similar to the (4+l)-dimensional Yang Mills model. 

The S 2 sigma model can also be regarded as the continuum limit of an 
array of Heisensberg ferromagnets. 

The static Lagrangian density for the S 2 sigma model is given by 

J |V ^' 2 ' 

where is a unit vector field. 

In the dynamic version of this problem, where <fi : IR 2+1 — ► S 2 , the La- 
grangian is 

L= I |^| 2 -|V0| 2 . 



Identifying S 2 = CP 1 = C U {oo} we can rewrite this as 

(4) 



\d t u\ 2 


\Vu\ 2 


(i + M 2 ) 2 


(i + M 2 ) 2 



The calculus of variations on this Lagrangian yields the following equation 
of motion for the S 2 sigma model: 

(1 + \u\ 2 )(d 2 u - d 2 x u - d 2 u) = 2u(\d t u\ 2 - \d x u\ 2 - \d y u\ 2 ). (5) 

Here u represents the complex conjugate of u. 

In this paper we first investigate the charge 1 sector of the (4+l)-dimensional 
Yang Mills model and then in the second part go on to the charge 2 sector 
of the S 2 sigma model. In each part, we first derive the evolution equation 
for our model. Second, we explain the numerical scheme used to investigate 
it. Third we go over the predictions made by the geodesic approximation for 
the model. Fourth, we go over the results generated from the computer runs. 
In both models the results encompass comparing the trajectory of the evolu- 
tion to that predicted by the geodesic approximation, and characterizing the 
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profile generated at fixed time slices. In both models our investigation of the 
profiles yields an improvement in approximate solutions of the differential 
equations, and these approximate solutions become exact in the adiabatic 
limit. 



2 The (4+1) -dimensional Yang Mills Model, 
Charge 1 Sector 

The first things to identify in this problem are the static solutions to equation 
(0). These are simply the 4 dimensional instantons investigated in |]J. In Jl| 
the form of all such instantons in the degree one sector is shown to be: 

, . , 1 f (x — a)dx — dx(x — a) 1 

A(x) = - < - — : \= } x = xi + x 2 i + x 3 j + x±k e H. 



2 [ A 2 + \x - a\ 
The curvature F of this potential is computed by F = dA + [A, A] and is: 

„ dx A dx 



(A 2 + |x-a| 2 ) 2 ' 



One notices the denominators are radially symmetric about x — a. 
An instanton of unit size centered at the origin would be 

A(x) 



2 [ l + \x\ 

One motion one can study from these static instantons would be to consider 
connections the form 

1 f xdx — dxx 



A{r,t) = 2 \ /(r,t)+r 2 J r = ^xl + x 2 2 + x 2 3 + x\ 

and derive an equation of motion for f(r,t) from (|3p. This is actually quite 
computationally extensive. 

To get at this, start with connections of the form 

A(r, t) = -g(r, t) {xdx — dxx} . 

Using the formula for F given by F = dA + [A, A], and a symbolic algebra 
program such as Maple, one can use (^j to calculate the differential equation 
for g: 
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%g = 12g 2 +^- + %g-8g 3 r\ (6) 
It is much less difficult to now compute a differential equation for f(r,t) 

if 

g\T,t) = . 

f(r,t)+r 2 

It is 

a?/ = #/ + - + (O,/) 2 - (a/) 2 ) ■ (7) 

The static solutions for f(r,t) are simply horizontal lines, f(r,t) = c. 
The geodesic approximation for motion under small velocities states that 
solutions should progress from line to line, i. e., f(r, t) = c(t). Here the length 
scale is given by \/c. f(r,t) = is a singularity of the system, where the 
instantons are not well defined. We observe progression from f(r, t) = Co > 
towards this singularity. The validity of the geodesic approximation can be 
evaluated from how f(0,t) evolves versus how it is predicted to evolve, and 
the differences between f(r,T) with T fixed and a horizontal line. 



2.1 Numerics for the (4+l)-dimensional Yang Mills 
Model 

In both the (4+l)-dimensional Yang Mills Model and the S 2 sigma model 
covered in section 3, a finite difference method is used to compute the evolu- 
tion of (0) and ( |12"D numerically. Unless otherwise noted, centered differences 
are used consistently, so that 

g(x + 5)- g(x - 5) 



d r g(x) 
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2 _ g(x + 5) +g(x-6) -2g(x) 

In order to avoid serious instabilities in (|7|), the terms 

dlf + ^ (8) 
r 

must be modeled in a special way. We see this revisited for the S* 2 sigma 
model in section 3. Allow 
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where 



£ = r ~ 5 d r r 5 d r . 



This operator has negative real spectrum, hence it is stable. However the 
naive central differencing scheme on (||) always results in uncontrolled growth 
near the origin. General wisdom holds that when one has difficulties with the 
numerics in one part of a problem one should find a differencing scheme for 
that specific part in the natural to that specific part. Applying this allowed 
for the removal of the problem near the origin. Instead of using centered 
differences on cftf and on d r f, we difference the operator: 



The "natural differencing scheme" is then 



f(r + S)-f(r) 
6 



r — 



S\ 5 ff (r)-f(r-6) 
6 



In |J, studying the charge 1 sector of the S 2 sigma model, stationary 
solutions were found to be unstable and to shrink spontaneously under the 
numerical scheme. This model exhibits no such difficulties. Unless bumped 
from the stationary state with an initial velocity, stationary solutions do not 
evolve in time. 

A discussion of the stability of this numerical scheme is available in [[3]. 

With the differencing explained, to derive f(r,t + At), one always has a 
guess for f(r,t + At) given by either the initial velocity, e.g. f(r,t + At) = 
f(r,t) + v At, or by f(r,t + At) = 2f(r,t) - f(r,t - At). Use this to 
compute d t f(r, t) on the right hand side of (^). Then solve for f(r,t + At) in 
the difference for d^f(r,t), and iterate this procedure to get a more precise 
answer. So, one iterates 



f(r,t + At) = 2f(r,t)-f(r,t-At) + (Aty 



d 2 r f(r,t) 



5d r f(r,t) 



2d t f(r,t) 2 2d r f(r,t) 2 8d r f(r,t)r 
f(r,t) + r 2 f(r,t) + r 2 f(r,t)+r 2 

where all derivatives on the right hand side are represented by the appropriate 
differences. 
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There remains the question of the boundary conditions. The function 
/ is only modeled out to a value r = R ^> 0. Initial data for /(r, 0) was 
originally a horizontal line. The corresponding boundary conditions are that 
f(R, t) = f(R - Ar, t), and that 

/(0,t) = |/(Ar,f)-±/(2Ar,t) 

i.e., that / is an even function. 

Subsequent investigation of the model indicated that the appropriate form 
for f(r, t) was a parabola instead of a line, and the f(R, t) boundary condition 
was changed to reflect this. For the runs with parabolic initial data, we set 
the boundary condition at R to be: 

d r f(R,t) = d r f(R-Ar,t) R 



R- Ar' 



2.2 Predictions of the Geodesic Approximation for the 
(4+l)-dimensional Yang Mills Model 

Equation (Q) gives us the Lagrangian for the general version of this problem. 
We are using 

A = - { — — — n — } r = \ x{ + X2 + XI + x%, 



2\ / + r 2 J 

and we use the geodesic approximation which predicts we move close to the 
moduli space of these solutions. Restricting the Lagrangian to this moduli 
space gives us an effective Lagrangian. The portion of the integral given by 



V^ijj Rij)dx 



represents the potential energy and integrates to a topological constant, hence 
it may be ignored. We see this same phenomenon again in section 3.2 in the 
study of the S* 2 sigma model. We need to calculate 



(d t Ai,dtAi)dx. 



First calculate 

(d t A i ,d t A i )= (/ + r2)4 - 
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So the Lagrangian is 

3r 2 (d t f) 2 



( / + r 2)4^ 

If we take / = constant, using y = x/y/f we can rewrite this integral as 



f j Ri (i + M 2 ) 4 



dy. 



The integral with respect to y converges, hence we have the effective La- 
grangian 

This is purely kinetic energy. Since the potential energy is constant, so is the 
kinetic energy. We have 

(W _ 

Integrating this we get 

/ = (cit + c 2 ) 2 . 
If / = occurs at time T, we find 

rp _ _ C 2 L 
Cl 

hence we rewrite this as 

f = a(t-T) 2 . 

This is how we predict that f(0,t) will evolve. This same evolution is pre- 
dicted in section 3.2. 



2.3 Results for the (4+l)-Dimensional Yang Mills Model, 
Evolution of f(0,t) 

The computer model was run under the condition that f(r, 0) = fo with 
various small velocities. The initial velocity is dtf(r, 0) = Vo, other input 
parameters are R = r max , Ar and At. 

The first question to ask is how does the evolution of the origin occur. We 
note that as r — > oo equation (|7D becomes the regular linear wave equation 

d 2 J = d*j, 
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and so the interesting nonlinear behavior is at the origin. Consequently we 
track the evolution of f(0,t). 

We find that the evolution of /(0, t) evolves exactly as predicted, as a 
parabola of the form: 

f(0,t) = a(t-T) 2 . 

Unsurprisingly, we obtain this same result in section 3.3 for the S 2 sigma 
model. 

The time to "blow up" is the parameter T in this equation. Recall that 
fo = f( r , 0) is the initial height and d t fo(r,0) = t> is the initial velocity. 
Using a least squares parabolic fit to the origin data obtained after f(0, t) < 
0.5fo, one obtains the parameters a and T for a given origin curve. Table |I| 
shows the behavior. 

We calculate the parameters a and T from the initial conditions, and we 
find 

v 2 
4/o' 

and 

j, _ 2/o 
\ v o[ 

A typical evolution of f{0,t) is given in Figure [I]. In this figure, the 
equation 0.000025(t — 200) 2 neatly overlays the graph of /(0, t). This picture 
represents the evolution where f = 1.0 and v = —0.01. Hence a = ( -°;° 1 2 = 
0.000025 and T = ^1 = 200, as predicted. 
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Table 1: Parabolic fit to /(0, t) vs. Initial conditions fo and vq 



fo 


vo 


a 


T 


1.0 


-0.010 


0.00002501 


200.1 


2.0 


-0.010 


0.00001257 


399.4 


0.5 


-0.010 


0.00005166 


99.0 


4.0 


-0.010 


0.00000626 


799.7 


4.0 


-0.020 


0.00002503 


400.1 


4.0 


-0.005 


0.00000157 


1599.3 
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2.4 Characterization of Time Slices f(r,T): Evolution 
of a Horizontal Line in the (4+l)-dimensional Yang 
Mills Model 

The most striking immediate result is that the initial line, /(r, 0) = /q, 
evolved an elliptical bump at the origin that grew as time passed. This 
happens again in section 3.4 with the S 2 sigma model. Figure |2| shows this 
behavior. 

The elliptical bumps can be modeled as 



a 



The question naturally arises as to how the parameters a, b and k evolve. 
This is straightforward: 

a = t 
k = fo + v t 
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Time Profiles of Data 




60 80 100 120 140 160 180 200 
r 



Figure 2: (4+l)-dimensional Yang Mills model, Time Slices f(r,T) evolve 
an elliptical bump at the origin 
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2.5 Characterization of Time Slices f(r,T): Evolution 
of a Parabola in the (4+l)-dimensional Yang Mills 
Model 

The elliptical bump that formed in the evolution of a horizontal line and the 
various configurations that ensued after it bounced off the r = and r = R 
boundary suggested that perhaps the curve was trying to obtain the shape 
of a parabola. After all, near r = 0, ellipses are excellent approximations for 
parabolas of the form 

f(r,t) = pr 2 + h. (10) 
To get the parabola, calculate from the general form of our ellipse in (|9|) 

dy x 2 b 2 
dx (y — k)a? 

so 

d 2 y —b 2 xb 2 dy 

dx (y — k)a 2 (y — k) 2 a 2 dx 
At x = 0, y — k = b and this gives 

d 2 y —b 
dx a 2 

Recall from the previous section that b = ct 2 and a — t, so this gives 

d 2 y 



The identification of c gives 



dx 2 



d 2 y _ v 2 



So 



dx 2 4/ 



1 d 2 y v 2 
P = -~-^= 



2 dx 2 8/ 

When a run is started with this initial data, d t fo = vo = —0.01, fo = 
/(0,0) = 1.0 and p = -gl = -0.0000125, the time slices of the data have 
this same profile. This is shown in figure H. 



13 



The curvature of the parabola at the origin, as measured by the parameter 
p from equation (|j~0D changes by less than 1 part in 100 during the course of 
this evolution, a graph of p over time can be seen in figure |j. 

The other parameter in equation (|10D for the parabola, h, should be given 
by the height of the origin, but this was calculated in the previous section to 
be a(t — T) 2 , substituting the expressions for c and T we obtain: 



h(t) = 4 f t - 2/ < 



4/o V M 

This is indeed the correct form, as shown in figure The initial conditions 
were v = -0.01 and f = /(0,0) = 1.0, hence h(t) = 0.000025(t - 200) 2 . 
The plot of the function overlays the data. 

Now, using both expressions for p and h, one can get the general form of 
a parabolic f(r,t), which is 

/M)= ^ + ii( t -^y (11) 

8/o 4/o V \vq\/ 

Substitute this into the partial differential equation (0), get a common de- 
nominator and simplify to obtain: 



4/o 



-^r 2 + 2 
4/o 



4/o 



2fo 
\ v o\ 



+ 2r 2 



4/o 



± r 2 + 2 ± 



4/o 



4/o 



2r 2 



Clearly these two sides are not the same, and the difference between them is 



4/o 

Since our concern is the adiabatic limit, Uq/(/o) i s always chosen to be small, 
and this difference goes to zero in the adiabatic limit. This term is always 
much smaller than the 2r 2 term, and since we only have accurate numerics 
away from the neighborhood of the singularity, the term with t — 2fo/\vo\ 
is on the order of magnitude fo/\vo\, which is large compared to the correc- 
tion. This provides an improved approximate solution for this model which 
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becomes exact in the adiabatic limit. This same behavior is found in section 
3.5 

The form of the connection A given by this parabolic model is 



A(r,t) = 




This is close to the connection A given by the geodesic approximation mul- 
tiplied by an overall factor. This form suggests a relativistic correction. 
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Time profiles with parabolic initial data 



t = 20s 
t = 40s 
t-60s 
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20 



30 



40 



50 



Figure 3: (4+l)-dimensiorial Yang Mills model: Time slices of the evolution 
of a parabola are parabolas. 
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Parabolic parameter p as a function of time 




40 50 60 
Time in seconds 



Figure 4: (4+l)-dimensiorial Yang Mills model: evolution of parabolic 
rameter p with time. 
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Evolution of parabolic parameter h with time 




20 40 60 80 100 120 140 160 180 200 

t 



Figure 5: (4+l)-dimensional Yang Mills model: Parabolic parameter h as a 
function of time evolves as f(0,t). 
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3 The Charge 2 Sector of the S 2 Sigma Model 



The first thing to identify in this probiem are the static solutions determined 
by equation (g). These are outlined in |J among others. The entire space of 
static solutions can be broken into finite dimensional manifolds M n consisting 
of the harmonic maps of degree n. If n is a positive integer, then M n consists 
of the set of all rational functions of z = x + iy of degree n. We restrict our 
attention to M2, the charge two sector, on which a generic solution has the 
form 

u = a + ; 

z z + oz + e 

depending on the five complex parameters a, f3, 7, 6, e. To simplify consider 
only solutions of the form 

1_ 

z 2 

with 7 real. The geodesic approximation predicts that solutions evolve close 
to 

tW 

In this investigation, we use a radially symmetric function f(r,t) = 7, and 
calculate the evolution equation for f(r,t). It is: 

a? / = d'tf + ^ - ^M. + JM_ ((dtff - (ojy) . ( i2) 

We can evaluate the validity of the geodesic approxmation by evaluating 
how f(0,t) evolves versus the predicted evolution and how f(r,T) with T 
fixed varies from being a horizontal line. 

Immediate similarities can be seen between ( |T2D and equation ([7p. The 
static solutions are f(r,T) = c for T fixed and any constant c. Here the 
length scale is given by yfc. We investigate progression from f(r, 0) = c > 
towards the singularity at f(r,T) = 0. 

3.1 Numerics for the S 2 Sigma Model Charge 2 Sector 

As with the other models, a finite difference method is used to compute 
the evolution of (|1^) numerically. Centered differences are used consistently 
except for 
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A naive central differencing scheme here yields serious instabilities at the 
origin. Similar to what we did in section 2.1, let 



-5;) ...-.;),. ( ;)2 t . 



Lf = r- b d r r b drf = d 2 J 

T 

This operator has negative real spectrum, hence it is stable. The natural 
differencing scheme for this operator is 



-5 



( f(r + S)-f(r) \ _ f_S\ 5 ff(r)-f(r-5) 



27 V 5 J V 27 V 5 



This is the differencing scheme used for these terms. 

In || studying the charge 1 sector of the S 2 sigma model it was found 
that under their numerical procedure stationary solutions were unstable and 
would evolve in time. This model has no such problems with the stationary 
solutions. They do not evolve in time unless first bumped with an initial 
velocity. 

Now with the differencing explained, we derive f(r, t + At) in exactly 
the same manner as for the (4+l)-dimensional model in section 2.1 and the 
charge 1 sector in ||]. We have an initial guess for f(r,t + At), either given by 
f(r,t + At) = f(r, t) + t>o At with vq the initial velocity given in the problem, 
or on subsequent time steps f(r,t + At) = 2f(r,t) — f(r,t — At). We use 
this guess to compute d t f(r,t) on the right hand side of (|T2]), and then we 
can solve for a new and improved f(r,t + At) on the left hand side of (|r2|) . 
Iterate this procedure to get increasingly accurate values of f(r,t). 

The boundary condition at the origin is found by requiring that f(r,t) is 
an even function, hence 

/(0,t) = ^/(Ar,t)-i/(2Ar,t). 

At the r = R boundary the function should be horizontal so f(R,t) = 
f(R — Ar, t). 

Subsequent investigation of this model indicated that the appropriate 
form of f{r,t) was a parabola instead of a line, and in order to investigate 
this phenomenon, the f(R, t) boundary condition was changed to reflect this. 



20 



For the runs with parabolic initial data, we set the boundary condition at R 
to be: 

d r f(R,t) = d r f(R-Ar,t)-^—. 

R — Ar 

3.2 Predictions of the Geodesic Approximation for the 
S 2 Sigma Model 

Equation (f|) gives us the Lagrangian for the general version of this problem. 
We are using 

A 

u = — 

z l 

for our evolution, and via the geodesic approximation we restrict the La- 
grangian integral to this space, to give an effective Lagrangian, as we did 
in section 2.2. The integral of the spatial derivatives of u gives a constant, 
and hence can be ignored. Under these assumptions, up to a multiplicative 
constant, the effective Lagrangian becomes 

r\d t \f 



L = rdr- ,, 
Jo (r 4 + A 2 ) 

which integrates to 

r . (Wtt 

Since the potential energy is constant, so is the kinetic energy, hence 

d t \ = kV\. 

Integrating this one obtains 

A = (cit + c 2 ) 2 . 
If A = occurs at time T, we find 

T - — — 

Cl 

hence we rewrite this as 

A(t) = a{t-Tf. 

This is exactly the same as the evolution predicted in section 2.2 for the 
Yang Mills Lagrangian. Since equation (|T2|) tends towards the linear wave 
equation when r —>■ oo, the interesting behavior occurs at the origin. This is 
how we predict that f{0,t) will evolve. 
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3.3 Results for the S 2 Sigma Model: Evolution of /(0, t) 



The computer model was run under the condition that f(r, 0) = fo with 
various small velocities. The initial velocity is dtf(r, 0) = v o, other input 
parameters are r max = R, Ar and At. 

A typical evolution of f(0,t) is given in Figure |6|. This is best modeled 
by a parabola of the form a(t — T) 2 , exactly as predicted. This curve in 
particular is approximated by 0.0000998(t — 100. ) 2 . 

We fit f{0,t) to a parabola of the form a(t — T) 2 for various initial con- 
ditions. Table |] gives initial conditions fo = f(r, 0) and v o = d t f(r, 0) and 
the subsequent T and a when Ar = 0.025 and At = 0.001. 

Exactly as in section 2.3, we have 




and 



4/o' 



Table 2: S 2 sigma model: Parameters for best fit parabola to f{0,t) vs. 
initial data fo and vq- 



fo 


^0 


T 


a 


1.0 


-0.01 


200 


0.0000250 


1.0 


-0.02 


100 


0.0000998 


1.0 


-0.03 


67 


0.000224 


1.0 


-0.04 


50 


0.000398 


0.5 


-0.01 


100 


0.0000499 


2.0 


-0.01 


400 


0.0000125 


3.0 


-0.01 


600 


0.00000833 
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Typical evolution of f(0,t) 




20 30 



40 50 60 70 

Time in seconds 



80 90 1 00 



Figure 6: S 2 sigma model: Evolution of /(0, t) and overlaying fit to parabola. 
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3.4 Characterization of Time Slices f(r,T): Evolution 
of a Horizontal Line in the S 2 Sigma Model 

With the evolution of /(0, t) taken care of, we consider the shape of the time 
slices f(r, T) for a given fixed T. As in section 2.4 with the (4+l)-dimensional 
Yang Mills model, this is rather striking. An elliptical bump forms at the 
origin, as seen in figure [7|. 



Time slices f(r,T) 




250 



Figure 7: S 2 sigma model: Time slices of f(r,T) evolve an elliptical bump 
at the origin. 



Once again exactly as in section 2.4 with the (4+l)-dimensional model, 
this elliptical bump has equation 

The question naturally arises how the parameters a, b, and k evolve, and as 
in section 2.4, we once again have 

a = t 
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k = fo + v t 



Close to the origin the ellipse is well approximated by a horizontal line, 
giving further credence to the validity of the geodesic approximation. 
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3.5 Characterization of Time Slices f(r,T): Evolution 
of a Parabola in the S 2 Sigma Model 

As with the (4+l)-dimensional model in section 2.5, the evolution of the 
ellipse suggested the curve was trying to obtain the shape of a parabola of 
the form 



f(r,t) = pr 2 + h. (14) 

To get the general form of the parabola, we follow the calculation from 
the (4+l)-dimensional model in section 2.5. From our ellipse equation dl3f) : 



2U2 



dy x b 



dx (y — k)a 2 

so 

d 2 y —b 2 xb 2 dy 

dx (y — k)a? (y — k) 2 a 2 dx 
At x = 0, y — k = b and this gives 

d 2 y -b 
dx a 2 

Recall from the previous section that b = ct 2 and a = t, so this gives 

d 2 y 



The identification of c gives 



So 

p = - 



dx 2 



d_y_ Vp 
dx 2 " 4/o' 

1 d 2 y v 2 



2 dx 2 8/ 

Rather unsurprisingly, this echoes the result of section 2.5. 

When a run is started with this initial data, d t fo — v o = —0.02, f = 

2 

/(0, 0) = 1.0 and p = = —0.00005, the time slices of the data have this 
same profile. This is shown in figure || The parabolic parameter p varies 
by less than 1 part in 10 during the run as seen in figure ^|. The parabolic 



parameter h evolves close to /(0, t) as seen in figure 10 
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For a parabola of the form 

f(r,t) = p(t)r 2 + h(t), 

we have 

v 2 

i.e. p(t) is constant, and 



4/o V M 

Using the identification of p and h in the parabolic form of f(r,t) and 
letting 

r = t _2/o 
\ v o\ 

we have 

f( r t ) = _il_ r 2 + ^o_ r 2 

Substitute this into the partial differential equation (|i~2|), get a common de- 
nominator, and simplify to obtain 

4 \ 2 

r 2 2 1 4 \ 1 V 4 JL 



32ff vJ -rr + rj + -r 



v o rA V 4 , V l 4 



32/ 3 4 + 32/ 3T + 2f r • 



The difference between the two sides is 



4 2 2 



64/ 3 32/ ( 







As in section 2.5 for the (4+l)-dimensional Yang Mills model, our concern 
is with the geodesic approximation, and so fg/(/o) is always chosen to be 
small. The correction is then much smaller than the term 



2/0 
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Figure 8: S 2 sigma model: Time slices of the evolution of a parabola are 
parabolas. 
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Parabolic parameter p as a function of time 




40 50 60 
Time in seconds 



100 



lire 9: S 2 sigma model: Evolution of parabolic parameter p with time. 
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Parabolic parameter h as a function of time 




-0.2 



10 20 30 40 50 60 70 80 90 100 

Time in seconds 



Figure 10: S* 2 sigma model: Evolution of parabolic parameter h with time, 
comparison to /(0, £). 
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4 Conclusions 



The formation of singularities in the geodesic approximation to the (4+1)- 
dimensional Yang Mills model and the (2+l)-dimensional S* 2 sigma model has 
been studied numerically using radially symmetric solutions and an iterated 
finite differencing scheme. 

The predictions of the Lagrangians of the two models are that the tra- 
jectory towards blow-up should occur parabolically as a(t — T) 2 where T 
represents the blow-up time. This is confirmed in the behavior of both nu- 
merical models. 

The geodesic approximation in both models predicts that the solution 
will evolve in time from horizontal line to horizontal line. To first approxi- 
mation, near the origin, this is what occurs. A more precise characterization 
is available. Both of these models, when started with a horizontal line as an 
initial condition, evolve an elliptical bump at the origin. 

The elliptical bumps suggested that the model preferred a parabolic initial 
condition and a parabolic state for the fixed time profile of the evolution. 

These parabolic solutions become exact in the adiabatic limit, and provide 
alternative approximate solutions to the differential equations. 
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